Circular hydraulic jump in generalized-Newtonian fluids 
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We carry out an analytical study of laminar circular hydraulic jumps, in generalized-Newtonian 
fluids obeying the two-parametric power-law model of Ostwald-de Waele. Under the boundary-layer 
approximation we obtained exact expressions determining the flow and an implicit relation for the 
jump radius is derived. Corresponding results for Newtonian fluids can be retrieved as a limiting 
case for the flow behavior index n = 1, predictions are made for fluids deviating from Newtonian 
behavior. 
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When a smooth jet of liquid falls vertically on a rigid 
horizontal boundary, it spreads out radially in a thin 
layer until a sudden increase in depth (of a circular ge- 
ometry) may occur; this is a Circular Hydraulic Jump 
(CHJ) (see Fig.l). CHJ is an observable in everyday 
life as one opens a water tap into a kitchen sink. Apart 
from this simple day to day observation, striking analo- 
gies with CJHs have been seen at femto-scale where a 
CHJ appear as metal femtoliter cups [ij, and also at as- 
tronomical scales where it is manifested as a white hole 
horizon [2]. 

In this letter we present new results on CHJs in 
generalized-Newtonian fluids, where, the effective viscos- 
ity depends, in general, on shear rates of the flow, which 
leads us to interesting results applicable to a larger class 
of fluids, results for Newtonian fluids can be obtained as 
a special case. In contrast to our study, earlier studies 
on CHJs has been done only for the Newtonian fluids. 

CHJs were noticed and theoretically investigated for 
the first time by Rayleigh in 1914, who developed a the- 
ory consistent for inviscid flow by conserving mass flux 
and momentum flux across the jump However, later 
studies revealed that due to divergence in the flow, the 
flow layer becomes very thin with increase in radial dis- 
tance, where viscosityplays a significant role in the flow 
dynamics 0, d, @, 0, Watson in 1964 considered the 
influence of fluid viscosity on CHJs Q ■ In the thin layer 




of the flow, boundary-layer theory was applied to model 
the flow 0. At the jump, application of the momentum 
theorem gave a prediction for the jump radius. Experi- 
ments of Watson Q , Olson and Turkdogan [5|] , Craik et 
al. [6], and Liu and Lienhard [lo| characterized the do- 
main in which Watson's theory showed good agreement 
with experiments, at the same time the region where the 
theory seemed weaker was also marked. The agreement 
was good for smaller jump height to jump radius ratio 
and poor in the opposite limit of relatively, large jump 
heights or small jump radii. Origin of the discrepancies 
was found in Watson's simplifying assumptions: unidi- 
rectional flow beyond the jump and neglecting surface 
tension effect. 

Several experimental and theoretical studies since then 
have revealed that CHJs have broadly two types of steady 



FIG. 1: A circular hydraulic jump (CHJ) is formed when a jet 
of liquid falls vertically on a rigid horizontal boundary (Photo 
courtesy of John Bush and Jeff Aristoff, MIT) 



states, classified as type I and type II [III, ll2|, |13| . In ex- 
periments, the depth H beyond the jump is controlled 
by a circular boundary of variable height H, sufficiently 
far from the jump. Type I steady state is observed below 
a certain critical value He, beyond which transition to 
type II steady state occurs. In type I steady state, where 
surface flow is unidirectional and boundary-layer sepa- 
ration occurs beyond the jump, Watson's assumptions 
were found to be appropriate. On the other hand type 
II steady state is distinguished by reverse surface flow 
(roller) adjoining the jump, the boundary layer separa- 
tion is pronounced at the jump as a result of which the 
flow need not decelerate significantly as it passes through 
the jump: the discontinuity in the radial speed assumed 
in Watson's work need not arise. 

Steady jumps of type II is an active area of research 
tiU date, Bohr et al [Jlil, Higuera [11], Yokoi [12] stud- 
ied the boundary layer separation beyond the jump, and 
developed explicit scaling law for the jump radius which 
conforms well with the experiments. The key feature 
of these studies has been analysis of the complex flow 
dynamics near the jump, and its effect on the jump ra- 
dius. In contrast to above approaches, Watson's theory 
gives an implicit relation for the jump radius consider- 
ing a simplified average behavior of the flow just after 
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the jump. Naturally as a consequence of this simplifica- 
tion, Watson's theory has limitations, typically when the 
boundary-layer separation is dominant near the jump viz. 
steady type II jumps. However, experiments suggest that 
in all types of CHJs: type I; type II or unsteady, Wat- 
son's prediction for the jump radius provides an adequate 
leading-order description. In the framework of Watson's 
theory, in 2003 Bush and Aristoff considered the surface 
tension effect in Type I jumps, which led them to a re- 
vised estimate for the jump radius [isj . It was shown 
that in type I steady jumps of small radii and heights, 
surface tension plays a significant role in the dynamical 
balance at the jump. An expression was derived for the 
radial force per unit length associated with the curvature 
at the jump and Watson's result for the jump radius was 
appropriately revised. Experiments done by the same au- 
thors |13|] showed good overall agreement with the theory 
for type I steady jumps, in particular experimental data 
for jumps of small radius and height, which were poorly 
described by Watson's result, matched well on including 
the surface tension correction. 

CHJs so far has been studied (both theoretically and 
experimentally), only for the Newtonian fluids, charac- 
terized by a viscosity coefficient, which does not depend 
on the shear rate of the flow. Many fluids in nature 
e.g., oil, honey, paints, polymeric solutions, lime-water 
mixture etc, are non-Newtonian, having an effective vis- 
cosity, which depends on the shear rates of the flows. 
Boundary-layer theory in non-Newtonian fluids has been 
extensively studied by James and Dabrowaski, Filipussi 
et al. UM- In this paper, we have modeled and stud- 
ied the CHJ in generalized-Newtonian fluids (a class of 
non-Newtonian fluids) obeying two-parameter power-law 
model of Ostwald-de Waele [l3|- In this model, for 
one dimensional flows, shear stress, t can be given by 
T = X (1^) , where parameters K and n are the flow 
consistency index and flow behavior index respectively 
(0 < n < 1 for a shear-thinning or pseudo-plastic fluid, 
for Newtonian fluids n — 1 and n > 1 for a shear- 
thickening or dilatant fluid). ^ is the shear rate or 
the velocity gradient perpendicular to the plane of shear. 
The quantity Veff — ^ (f^)" represents an apparent 
or effective viscosity as a function of the shear rate. Con- 
trary to the inflnite range of effective viscosity suggested 
by the mathematical form of the power-law model, in 
real fluids effective viscosity is bounded for all possible 
shear rates. Therefore, the power-law is approximate, 
and is only a good description of fluid behavior across 
the range of shear rates to which coefficients are fitted. 
Nevertheless power-law is the simplest and most widely 
used model to describe the non-Newtonian fluid behav- 
ior, permit mathematical predictions, and correlate ex- 
perimental data. 

Our study is in the framework of Watson's theory in 
a sense that, assuming a continuous variation in flow 
behavior from a Newtonian to non-Newtonian fluid, we 




FIG. 2: A schematic illustration of the boundary layer struc- 
ture establised within the CHJ. The viscous boundary layer 
grows until it reaches the surface at a radial distance r;,, from 
the point of impact, rj is the jump radius. 



have derived our results based on similar approximations 
as that of Watson Q . Laminar flow in thin layer is stud- 
ied under the boundary-layer approximation. Descrip- 
tion of the flow is given in terms of Blasius sub layer 
developing near the point of impact and a far-fleld simi- 
larity solution. The similarity-profile of the far field sim- 
ilarity solution has been applied in the Blasius sub layer, 
along with the standard Karman-Pohlhausen approxima- 
tion f^. Momentum theorem is applied including the 
surface tension correction term introduced by Bush and 
Aristoff [l3|, to obtain a prediction for radius of jump in 
power-law fluids. 

On considering a smooth jet of fluid descending verti- 
cally on a horizontal plane covered by the same fluid of 
depth H; from the point of impact the fluid spreads radi- 
ally in a thin layer, before the jump occurs. A boundary 
layer grows from the stagnation point up to a radial dis- 
tance ri, beyond which, it absorbs the whole of the flow. 
The jump occurs at some radial distance Vj. Let, a be 
the radius and Uq be the velocity of the jet ( just before 
the impact) , Q be the volume flux of the flow, U (r) the 
free surface stream velocity, and h(r) height of free sur- 
face from the bottom, r being the radial distance from 
the point of impact (Fig. 2). For r < r^,, the flow speed 
outside the boundary-layer ( of thickness S(r) ) remains 
almost constant, equal to Uq, as the fluid here is unaf- 
fected by the viscous stresses. For r > n, the viscous 
stresses become appreciable right up to the free surface, 
the whole flow is of the boundary-layer type. 

For a stationary, radially symmetric flow with a free 
surface, boundary-layer equations for a power-law fluid 
are: 
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where r, z are cylindrical co-ordinates, with z mea- 
sured vertically upwards from the plate, and u(r, z), 
w{r, z) are the corresponding velocity components. 
Gravitational pressure gradient has been neglected 
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as in The boundary conditions are no-slip on 

the bottom: u(r, 0) = w(r, 0) = 0, no stress on the 
free surface: ^ = 0, and the kinematic boundary 
condition at the free surface: w{h,r) — u{h,r)^, 
which gives condition of constant volume flux 
Q = 2Tir udz = TTo^Uo =constant. Due to sym- 
metry of the flow shear stress components t^^ and r^g 



are zero, shear stress component Tro — K 



a{u) 



-d{u) 



In the region r > ri,, assuming far field solution of a 
similarity form: u = U{r)J{rf) where -q — , boundary 
conditions transform to: /(O) = 0, /(I) = 1, /'(I) = 0, 
and condition of constant volume flux can be written 
as: Q — 2'iTrUh f{ri)dr] = constant. On solving the 
equations of motion ([T]) and ^ for similarity solution, 
we found that similarity function is given by the implicit 
relation 

(3) 
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where the constant: c = |/3(-|, Here /? is the beta 

function and f{il)d7] evaluates to 3j/3(|, (say, = 
A), free surface velocity U{r) and free surface height h{r) 
for r > rii, turns out to be: 
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the constant, fc — 3 ifc^"^^^ and I is a length 

scale to be determined. For the determination of length 
scale we proceed by considering the flow in the region 
r < rb- The boundary-layer equation can be further sim- 
plified by using a standard averaging technique of von 
Karman and Pohlhausen Q , we shall use the self-similar 
velocity profile: u — C/(r)/(|), here 6 is the boundary- 
layer thickness. As a result of averaging over z, the mo- 
mentum integral equation of the flow is: 




{Uqu — v?)dz 




(6) 



solving it for the self-similar velocity profile leads to, 
'B{n + l) 



Sir) 
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where the constant B ~ — - — zr^pr 
the condition of constant volume flux at r 
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From here applying 
n gives. 



(8) 



matching U{r) = Uq a.t r = ri, leads to: 

' [(1 - 2Ac)n + 1 + Ac]{n + 2)7r"~2 
(2n- l)(n-M)c"+i(2A)"+i 

where, Re — j^^i-in , is the jet Reynold's number. 
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FIG. 3: Variation of boundary layer thickness 5, free surface 
height h and free surface speed U with radial distance, for 
the flow behavior indices n = 0.8, 1.0, 1.2. Propagation of 
viscous effect to the free surface, free surface height decay 
and momentum loss in the flow is delayed with increase in 
flow behavior index n. Results of Watson are retrieved for 
n = 1. 
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On applying constant volume flux condition Q = 

, free surface height in 
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the region r < can be obtained as: 
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The above expressions (jH), d?]), (fTO| . for free sur- 
face stream velocity U{r), free surface height h{r) and the 
boundary-layer thickness 6{r) in non-Newtonian power- 
law fluids, which leads to the prediction of jump radius, 
are some significant results of our paper. Fig. 3 shows 



the variation of 



U{r) h(r) 



^i?e^ and ^i?e^, as func- 



tions of -i?e"+2, for flow behavior indices 0.8, 1.0, 1.2. 
It can be observed that, propagation of viscous effect to 
the free surface, free surface height decay and momen- 
tum loss in the flow is delayed with increase in flow be- 
havior index n. Due to the viscous effect fluid within the 
boundary-layer loses its momentum as it proceeds up- 
stream, so the shear-rate decreases with increasing r. In 
a shear-thinning fluid viz. n = 0.8 effective viscosity in- 
creases as the flow progresses, viscous effect reaches the 
free surface earlier. If experiments are done flxing the jet 
Reynolds number jump radius is expected to be smaller 
then in a Newtonian fluid [n — 1) . Opposite behavior is 
expected for shear-thickening fluid in comparison to the 
Newtonian counterpart. 

For deriving jump radius, Watson applied momen- 
tum theorem at the jump, by equating the difference 
in pressure force across the jump to the difference in 
radial momentum flux. Bush and Aristoff [l3| showed 
that for circular jumps a radial force due to the sur- 
face tension effects arises at the jump. They considered 
the geometry of the jump and derived a curvature force 
Fc = — 27r(TAi/, where a is coefficient of surface tension 
and Ai/ = H — h(rj). By including this correction into 
the momentum equation they modifled the result of Wat- 
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where Bq = ^^^^ — is the jump Bond number. 

As the jump may occur at any point in the develop- 
ment of the flow. For a power-law fluid we obtain that: 
If the jump occurs after ri, i.e., rj > ri,, 
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for the jumps occurring before rh i.e., rj < rh, 



(12) 



rh(rj) 
JO 



u^dz 



h{rj 



1 



3nc 



iM(^.) 



Ul (13) 



So we get an implicit algebraic expression for the jump 
radius rj on substituting equations p2|) and ([13]) in (fTTj) . 
Here the flow rate Q, jet radius a, jump height H, coeffi- 
cient of surface tension cr, flow consistency index K and 
flow behavior index n are known physical quantities. 



In summary, we have studied laminar CHJs in non- 
Newtonian power-law fluids. Assuming a continues vari- 
ation in the flow behavior from Newtonian to non- 
Newtonian fluids, we have generalized the results of Wat- 
son Q for Newtonian fluids, to generalized-Newtonian 
fluids obeying the power-law model. Expression for the 
jump radius is obtained, which gives new predictions for 
power-law fluids. Our result reduces to that of Watson 
for flow behavior index n = 1. Surface tension correc- 
tion term introduced by Bush and Aristoff [3], which 
becomes significant for jumps of small radius and height 
has been taken care of. 
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